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Abstract 

Mixing  of  primary  cold  flow  air  and  secondary  helium  to 
control  the  ratio  of  specific  heats  for  the  medium  flowing 
through  a  simulated  chemical  laser  nozzle/lasing  cavity  was 
accomplished.  The  effects  of  a  range  of  mixture  specific 
heat  ratios  on  flowfield  behavior  were  examined  using  static 
pressure  ports  in  the  test  cavity.  Schlieren  photography 
and  high  speed  filming  aided  description  of  the  flow 
dynamics.  Results  indicated  that  boundary  layer  effects 
became  evident  in  the  nozzles  as  specific  heat  ratios 
increased.  Large  pressure  fluctuations  were  observed  in  the 
cavity  when  helium  was  introduced  into  the  flow  to  raise  the 
specific  heat  ratio.  This  unstable  behavior  was  attributed 
to  the  helium  mass  flow  into  the  mixer  and  the  mixer  design 
itself.  Use  of  the  air/helium  mixer  brought  about  the 
pressure  fluctuations  earlier  in  a  test  run  than  with  the 
mixer  removed  under  the  same  conditions.  Favorable  pressure 
conditions  for  lasing  were  achieved  for  at  least  two  seconds 
for  the  supersonic  nozzles'  design  specific  heat  ratio  of 
1.51.  Adverse  pressure  behavior  was  also  attributed  to 
three  dimensional  viscous  effects  along  the  cavity  walls. 
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EFFECTS  OF  SPECIFIC  HEAT  RATIO  ON  A 


SIMULATED  CHEMICAL  LASER  CAVITY  FLOW 


L: _ Introduction 


Background 

Interest  in  chemical  laser  research  has  increased  due 

the  possible  application  of  this  type  of  laser  to  the 

Strategic  Defense  Initiative  (SDI) .  The  Alpha  chemical 

laser,  a  key  component  of  the  SDI  Zenith  Star  program,  is 

now  in  development  and  testing  (1:23).  The  decision  on 

whether  or  not  the  chemical  laser  will  be  utilized  is 

related  to  the  performance  of  the  laser  which  depends 

directly  on  the  flow  structure  in  the  lasing  cavity  (2:1-1). 

Patterson,  Batten,  and  Howie  introduce  their  source 

flow  chemical  laser  with  the  following: 

Figure  1  shows  a  typical  source  flow  laser 
configuration.  These  lasers  are  characterized  by 
supersonic  oxidizer  streams  which  enter  the  laser 
cavity  from  an  array  of  source  flow  nozzles.  A 
secondary  fuel  stream  is  injected  either  sonically 
or  supersonically  into  the  oxidizer  stream  through 
wedges  placed  at  the  exit  of  the  source  flow 
nozzle.  The  mixing  of  the  two  streams  and  the 
subsequent  reaction  generates  vibrationally 
excited  molecules  which  are  the  source  of  the 
laser  radiation.  As  indicated  by  the  figure,  the 
plane  in  which  mixing  occurs  is  orthogonal  to  the 
expansion  plane  and  lasing  [cavity]  axis.  The 
purpose  of  the  source  flow  is  to  reduce  the  amount 
of  diluent  required  by  using  area  expansion  to 
help  control  the  temperature  and  pressure  rise  in 
the  laser  cavity.  (3:1) 
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Figure  l.l  Source  Flow  Chemical  Laser  Concept  (3:1) 


Laser  performance  is  strongly  connected  to  the  fluid 
mechanics  processes  of  the  flow  crossed  by  the  beam. 
Rapagnani  and  Lankford  from  the  Air  Force  Weapons  Laboratory 
(AFWL)  state 


Computer  modelling  of  the  flow  through  cavities  is 
not  a  trivial  task.  Analysis  of  the  flows  are 
complicated  by  the  existence  of  strong  cross 
stream  pressure  gradients  and  subsonic 
recirculation  zones  that  are  embedded  in  the 
supersonic  cavity.  This  analysis  of  a  high  base 
relief  cavity  has  indicated  the  existence  of  large 
property  variations  and  significant  regions  of 
sxibsonic  flow.  These  subsonic  recirculating 
regions  are  essentially  dead  spaces  in  the  cavity 
since  mass  is  neither  convected  in  or  out  of  an 
area  surrounded  by  a  closed  streamline.  These 
regions  are  local  hot  spots  since  the  temperature 
is  approaching  its  stagnation  value.  Large 
pressure  variations  occurred  in  this  cavity. 
(3:344) 
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The  previous  discussion  illustrates  the  great  need  for 
experimental  investigation  of  chemical  laser  cavity  flow- 
fields.  Research  is  needed  not  only  to  better  "map  out"  the 
quantitative  and  qualitative  structure  of  the  flow  in  the 
lasing  cavity;  but,  also  to  make  available  data  which  can  be 
used  for  comparison  of  results  obtained  through  computer 
simulation. 

Scope 

The  thesis  research  completed  by  Stiglich  (4)  is  the 
basis  of  this  investigation.  A  major  portion  of  his  work 
involved  the  design  and  construction  of  a  test  facility 
capable  of  investigating  flow  conditions  in  a  simulated 
laser  nozzle  assembly  and  lasing  cavity.  The  nozzle 
assembly  was  modelled  after  the  Alpha  laser  which  was 
designed  for  use  with  a  gas  specific  heat  ratio,  y,  of  1.51. 
Stiglich 's  study  used  only  air  as  the  medium  which  has  a  y 
value  of  about  1.4.  Using  a  cold  flow  mixture  of  compressed 
air  and  bottled  helium,  a  method  of  controlling  the  gas  y 
value  was  devised  for  this  investigation.  The  addition  of 
hypersonic  wedges  (see  Figure  1.1)  for  secondary  flow 
injection  was  not  accomplished  in  this  research.  Based  on 
time  histories  of  pressures  in  the  cavity,  schlieren 
photographs  of  the  cavity  flowfleld,  and  high  speed  film  of 
the  flow  dynamics,  the  operating  characteristics  of  the 
nozzle/cavity  system  of  this  simulated  section  of  a  chemical 


1-3 


laser  were  investigated.  The  data  was  scrutinized  for 
patterns  and  effects  caused  by  use  of  different  y  values. 

Additionally,  the  Data  Acquisition  System  (DAS)  used 
for  this  research  was  substantially  upgraded  to  include  20 
channel  simultaneous  recording  of  pressure  transducer 
voltage  output  throughout  a  single  test  run.  Previous 
equipment  allowed  only  6  pressure  readings  to  be  evaluated 
per  run.  Data  reduction  software  was  developed  to  convert 
transducer  voltage  output  to  pressure  values,  calculate  mass 
flow  rates  of  the  primary  and  secondary  gases,  determine 
mixture  y  values,  and  compute  nozzle  exit  Mach  numbers. 
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11^ _ Theory 


Analysis  of  the  pressure  data  obtained  during 
experimentation  included  determination  of  the  mixture 
specific  heat  ratio  for  the  combination  of  air  and  helium, 
and  nozzle  fluid  dynamics. 

Mixture  Specific  Heat  Ratio 

Several  assumptions  were  made  to  calculate  the  ratio  of 
specific  heats,  y,  for  the  gas  after  the  air  and  helium 
supply  mixed.  Since  cold  flow  was  used,  the  inherent 
complexity  of  high  temperature  reaction  chemistry 
common  in  actual  chemical  laser  operation  was  greatly 
simplified.  The  gases,  air  and  helium  at  room  temperature, 
were  assumed  to  behave  as  perfect  gases  and  applicable 
-elations  were  employed.  The  specific  heats  of  the 
individual  supply  gases  were  considered  constant  in  the 
temperature  range  of  testing.  Wark  (5:844)  states  that  the 
values  of  specific  heats  for  monatomic  gases  over  a  wide 
range  of  temperatures  are  constant.  Table  2.1  lists  the 
specific  heat  values  used  for  air  and  helium. 


Table  2.1  Specific  Heat  Values 


Gas 

(Molecular  Wt.) 

S 

(Btu/lb^-'F) 

Cy 

(Btu/lb||^-F) 

Air  (28.97) 

0.240 

0.17155 

Helium  (4.003) 

1.24906 

0.74944 
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The  ratio  of  specific  heats  is  defined  as 


To  determine  y  for  the  mixture  of  gases  the  conservation  of 
mass  yields 


•  =  •  +  • 


For  the  mixture  specific  heats 


•  ^  • 


Combining  the  above  equations  yields 


^Pmix  _  ^aii^Pair  ^Ha^PHe 


2 
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So,  given  the  mass  flow  rate  of  each  gas,  the  mixture  y  can 
be  determined.  The  mass  flow  rate  for  the  air  line  was 
determined  using  a  square-edge  orifice  meter  and  for  the 
helium  line  via  a  venturi  flow  meter.  The  calculation 
methods  for  these  meters  are  described  in  reference  6. 


Nozzle  Flow 

Since  the  converging-diverging  laser  nozzles  were 
contoured,  transverse  velocities  inside  were  considered  to 
be  small  and  the  flow  was  treated  as  one-dimensional.  From 
Shapiro  (7:83),  the  following  relation  is  given  for 
isentropic  flow  of  a  perfect  gas: 


£o 

p 


(6) 


Solving  for  Mach  number  and  applying  to  the  nozzle  exit 
conditions  gives  an  equation  for  the  nozzle  exit  Mach  number 
based  on  experimental  data 


2 

Y-1 


(7) 


where  p^  is  the  stagnation  pressure  upstream  of  the  nozzles, 
p^  the  static  pressure  at  the  nozzle  exit  and  y  the  value 
obtained  from  equation  5.  Equations  5  and  7  were 
incorporated  into  the  data  reduction  computer  program 
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POSTRUN. BAS  listed  in  Appendix  B. 

A  theoretical  counterpart  for  the  exit  Mach  number  of 
the  nozzles  was  derived  from  the  isentropic  equation  given 
in  Shapiro  (7:86)  for  area  ratio  as  a  function  of  y  and  Mach 
number  applied  to  the  exit 


=  c  = 


(T-l) 


(8) 


Solving  for  M^  gives  a  theoretical  equation  for  Mach  number 
as  a  function  of  y  and  the  constant  area  ratio  for  the  laser 
nozzle  of  32.6 


*2  (9 

M  =  e)  -  2 

*  Y-l 

Solving  M^  iteratively,  theoretical  and  experimental  exit 
Mach  numbers  for  a  range  of  y  values  were  then  plotted  for 
comparison. 


2 


4 


Ill .  Bxparinental  Apparatus 


Four  major  systems  make  up  the  experimental  apparatus: 
the  flow  system,  the  test  section  nozzle/cavity  system,  the 
data  acquisition  system,  and  the  schlieren  optical  system. 

Flow  System 

Figure  3 . 1  diagrams  the  flow  system  for  the  blowdown 
wind  tunnel.  Description  of  each  of  the  system  components 
follows: 

Compressed  Air  System.  The  compressed  air  system 
consists  of  two  ATLAS  COPCO  Air  Compressors,  Model  GAU-807; 
two  DELTECH  Filters,  Model  819;  two  Pioneer  Refrigerant  Air 
Dryers,  Model  R550A;  and  two  Arrow  Compressed  Air  Dryers. 
Each  of  the  air  cooled  compressors  is  capable  of  providing 
520  standard  cubic  feet  per  minute  (SCFM)  at  102  pounds  per 
square  inch  gauge  (psig) .  The  DELTECH  filters  are  used  to 
eliminate  any  pollutants,  such  as  oil,  from  the  air  flow. 
The  Pioneer  refrigerant  air  dryers  eliminate  moisture  from 
the  air  flow  and  support  a  mass  flow  rate  up  to  500  SCFM. 
The  calculated  pressure  loss  across  the  refrigerant  air 
dryers  is  approximately  3.6  psig.  A  second  drying  stage  is 
available  through  the  Arrow  Chemical  Air  Dryers.  Normally, 
the  chemical  air  dryers  are  by-passed  due  to  their  large 
pressure  loss,  approximately  10  psig.  Once  the  air  flow  is 
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through  the  final  stages  of  drying,  it  is  piped  to  the 
laboratories.  In  the  laboratory  the  air  goes  through  a 
globe  valve  and  into  a  cyclone  separator.  The  cyclone 
separator  is  the  final  station  used  to  eliminate  particulate 
matter  from  the  air  flow.  Next,  the  air  is  routed  through 
an  electrically  controlled  3  inch  ball  valve.  The  Starlite 
model  666-F  Ball  Valve  provides  quick  opening  via  pneumatic 
actuation  and  the  3  inch  diameter  provides  the  maximum  air 
flow  to  fill  the  horizontal  stilling  chamber  and  start  the 
nozzles  (4 : 3-1) . 

Orifice  Meter.  Mass  flow  measurement  of  the  compressed 
air  was  calculated  using  a  scpjare  edge  orifice  meter 
with  a  0.8  inch  throat.  Calculation  of  the  mass  flow  rate 
was  accomplished  using  the  method  of  reference  6. 

Horizontal  Stilling  The  six  foot  horizontal 

stilling  chamber  houses  a  steel  honeycomb  mesh  to  help 
straighten  and  distribute  the  air  evenly  prior  to  entering 
the  mixer. 

Helium  Supply  System.  The  supply  of  secondary  helium 
was  provided  from  a  farm  of  five  helium  tanks  at  2000  psig 
pressure.  Pressure  reduction  was  provided  by  two  Grove 
Powreactor  Dome  Controllers.  The  first  dome  valve  was  a 
Model  WBX-304-K3  and  reduced  the  storage  pressure  of  2000 
psig  to  500  psig.  The  helium  then  passed  through  a  venturi 
flow  meter  followed  by  the  second  dome  valve,  a  Model 
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Figure  3.1  Flow  System  Layout 


11010-P055B,  which  further  reduced  the  flow  pressure  to  the 
desired  pressure.  Exiting  this  dome  valve,  the  flow 
encountered  a  solenoid  valve  which  was  controlled  from  the 
operator  control  panel.  With  this  valve  open,  the  helium 
flow  is  choked  via  a  control  nozzle  which  was  designed  for 
the  purpose  of  preventing  pressure  fluctuations  downstream 
of  the  nozzle  from  varying  the  helium  mass  flow  rate. 

Figure  3.2  depicts  the  design  of  this  nozzle  which  consisted 
of  a  brass  insert  silver  soldered  into  a  19/32  inch  inner 
diameter  stainless  steel  pipe  connector.  The  brass  was  then 
drilled  to  the  desired  throat  diameter  and  the  nozzle 
shaped.  Once  through  this  flow  controller,  the  helium 
entered  the  horizontal  stilling  chamber  and  the  mixer. 
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Venturi  Flow  Meter.  Mass  flow  measurement  of  the 


helium  supply  was  accomplished  using  a  venturi  flow  meter 
with  a  0.5  inch  diameter  throat.  The  mass  flow  rate  was 
calculated  using  the  method  of  reference  6. 

Mixer.  The  mixer  design  was  based  on  the  Vassilatos 
parallel  multitube  ejector  design  described  by  Zakanycz  (8) . 
The  dimensions  of  the  discharge  end  of  the  Vassilatos 
multitube  ejector  were  enlarged  so  that  a  large  scale 
version  could  be  placed  inside  the  horizontal  stilling 
chamber.  The  mixer  was  made  of  stainless  steel  tubing, 
plywood,  Plexiglas,  8  inch  PVC,  and  aluminum  plate.  The 
mixer  design  is  shown  in  Figure  3.3.  The  mixer  slipped  into 
the  forward  section  of  the  horizontal  stilling  chamber  and 
was  held  in  place  by  sandwiching  the  flange  end  of  the  mixer 
plate  between  the  flange  ends  of  the  forward  and  rearward 
portions  of  the  horizontal  stilling  chamber.  Entering  the 
mixer,  the  air  was  forced  to  flow  around  a  chamber  full  of 
glass  beads  (marbles) .  The  glass  beads  caused  an  even 
distribution  of  air  to  flow  into  the  air  tubes.  A  fine  wire 
mesh  v'3s  placed  between  the  glass  beads  and  the  entrance  to 
the  tubes  to  filter  out  any  debris.  Once  in  the  tubes,  the 
air  was  directed  through  the  remaining  part  of  the  mixer  and 
out  the  end  plate.  The  secondary  flow  of  helium  entered  the 
sidewall  of  the  horizontal  stilling  chamber  and  was 
deflected  by  the  sidewall  of  the  mixer  to  flow  into  holes 
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Figure  3.3  Air/Helium  Mixer  Design 


drilled  in  the  side  of  the  PVC  cylinder.  Once  inside  this 
cavity  in  the  mixer,  the  helium  flow  was  forced  to  move 
around  the  air  tubes  traversing  the  cavity,  into  a  second 
set  of  tubes,  and  then  directed  out  the  end  plate.  The 
alternating  of  the  primary/secondary  tubes  exiting  the  end 
plate  induced  a  thorough  mixing  in  the  downstream  end  of  the 
horizontal  stilling  chamber  (4:3-4). 

Variable  Supersonic  Diffuser.  A  supersonic  diffuser 
was  built  to  maintain  the  low  pressures  in  the  cavity 
portion  of  the  test  section.  The  diffuser  was  made  out  of 
aluminum  and  plexiglas  and  actuated  by  a  24  v  dc  Delco 
motor.  Power  for  the  motor  was  provided  by  the  laboratory 
dc  power  supply,  a  Rapid  Electric  Co.  Model  S-528.  The 
motor  was  electrically  controlled  from  the  operator  control 
panel.  Figure  3.4  shows  the  motor  mounted  on  the  outside  of 
the  diffuser  box.  This  motor  drove  a  shaft  that  penetrates 
the  diffuser  side  wall  through  a  pressure  tight  seal.  The 
shaft  drives  a  gearing  mechanism  that  turns  a  threaded 
shaft.  The  threaded  shaft  moves  back  and  forth  either 
pushing  or  pulling  the  diffuser  slide  bracket  which  in  turn 
closes  or  opens  the  diffuser  blades  (diminishing  or 
enlarging  its  throat) .  To  protect  the  apparatus  from 
damage,  contacts  positioned  on  the  threaded  shaft  deactivate 
the  motor  as  contact  bars  close  micro-switches  mounted  on 
the  adjustment  bracket  behind  the  gear  box  (4:3-6). 
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Figure  3,4  Supersonic  Diffuser,  Top  View 

Vacuum  System.  To  establish  the  desired  test  section 
operating  conditions  it  was  necessary  to  have  high  pressure 
on  one  end  of  the  wind  tunnel  and  a  vacuum  on  the  other  end. 
The  vacuum  system  consists  of  16  large  tanks  suspended  from 
the  laboratory  ceiling  and  three  vacuum  pumps  located  in  the 
pump  room.  Two  of  the  three  vacuum  pumps  are  Leiman,  no 
model  number,  radiator  cooled  pumps,  and  the  third  pump  is  a 
Stokes,  Model  212-H.  The  Stokes  piunp  has  the  capacity  to 
pump  140  SCFM.  The  16  vacuum  tanks,  along  with  the  piping, 
provide  approximately  650  ft^  of  vacuum  space.  A  network  of 
8  inch  PVC  pipe  connected  the  test  section  to  the  vacuum 
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system.  The  interface  between  the  PVC  pipe  and  the  test 
section  was  a  Consolidated  Vacuum  Corporation,  Model  VCS-41, 
high  vacuum  valve.  The  high  vacuum  valve  allowed  the 
operator  to  isolate  the  wind  tunnel  from  the  vacuum  system. 

A  Meriam  U-tube  mercury  manometer.  Model  20DA40,  was  incorp¬ 
orated  to  provide  the  operator  with  the  vacuum  pressure 
reading  while  the  tunnel  was  being  evacuated  (4:3-7). 

Nozzle/Cavitv  Test  Section.  The  nozzle/cavity  design 
selected  and  its  general  operating  conditions  were 
established  from  information  provided  by  the  Air  Force 
Weapons  Laboratory  (AFWL) ,  Kirtland  AFB,  NM.  Figure  3.5 
shows  the  design  of  one  nozzle  segment  and  Table  3.1  (9) 
lists  the  coordinates  from  the  throat  (0.0000,0.0000)  to  the 
exit  (0.5560,1.1530)  of  the  nozzle  contour.  The  nozzle 
design  parameters  are  listed  in  Table  3.2.  Three  pressure 
ports  were  drilled  into  each  of  the  nozzle  segments  to  allow 
static  pressure  measurements  in  the  region  of  the  base  and 
the  nozzle  exit  plane.  Figure  3.6  displays  the  three  nozzle 
segments  as  they  are  mounted  in  the  test  section.  The 
current  design  allows  enough  material  to  machine  the  base 
region  to  accommodate  future  modifications  for  attachment  of 
the  secondary  flow  wedges  which  were  not  included  in  this 
investigation  (4:3-9).  The  test  section  was  constructed  out 
of  aluminum  and  the  windows  were  cut  from  1  inch  thick 
optical  plexiglas. 
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Table  3.1  Nozzle  Contour  Coordinates 


Z  Position  (inches) 

Y  Position  (inches) 

1 

0.0000 

0.0000 

2 

0.0013 

0.0240 

3 

0.0055 

0.0480 

4 

0.0204 

0.0960 

5 

0.1429 

0.2880 

6 

0.2604 

0.5040 

7 

0.3657 

0.7200 

8 

0.4742 

0.9600 

9 

0.5560 

1.1530 

3-10 


Table  3.2  Nozzle  Design  Parameters 


Parameter 

Quantity 

Specific  heat  ratio,  y 

Upstream  stagnation  pressure,  P^j 

Upstream  stagnation  temperature, 

Isentropic  exit  Mach  number, 

Isentropic  exit  pressure, 

Nozzle  area  ratio,  A^A* 

1.51 

100  psia 

73®F 

6.21 

0.0886  psia 

32.589 

Figure  3.6  Test  Section  Nozzle  Configuration 


Data  Accniisition  System 


Figure  3.7  depicts  the  Data  Acquisition  System  (DAS) 
hardware.  Analog  signals  from  piezoresistive  pressure 
transducers  are  passed  to  Signal  Conditioning  Units  (SCUs) 
for  amplification.  From  the  SCU  each  analog  signal  is 
transferred  via  a  BNC  cable  to  the  Nicolet  System  500  Data 
Acquisition  System  which  consists  of  the  Data  Acquisition 
Unit  (DAU) ,  the  Controller,  and  the  Video  Display. 

Piezoresistive  Pressure  Transducers.  Nineteen 
piezoresistive  pressure  transducers  were  used  to  monitor  the 
flow.  Table  3.3  shows  the  position  number,  appointed 
Nicolet  channel,  location,  serial  number  (ENDEVCO  unless 
otherwise  stated) ,  linear  operating  range,  and  signal 
conditioner  assigned  to  each  transducer.  Figure  3.8 
displays  the  available  transducer  tap  locations  in  the  test 
section.  ENDEVCO  Models  8506-2,  8506B-5,  and  8510B-5  were 
used  in  the  test  section.  Their  low  range  of  linear 
operation  (the  dash  numbers  indicate  the  psig  range)  was 
needed  for  the  very  low  pressure  of  the  flow  beyond  the 
nozzle  exit.  An  ENDEVCO  Model  8510B-15  was  utilized  in  the 
vertical  stilling  chamber  to  monitor  the  evacuation  of  the 
system  prior  to  a  test  run.  An  ENDEVCO  Model  8510B-100  was 
placed  in  the  horizontal  stilling  chamber  to  measure  the 
stagnation  pressure  just  upstream  of  the  test  section 
entrance.  Pressure  measurements  were  also  needed  for  air 
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Figure  3.7  Data  Acquisition  System 


T^□>le  3.3  Transducer/Thenmocouple  Locations 


Position 

Number 

Nicolet 

Channel 

Location 

S/N  (Model) 

Range 

(psig) 

SCU  No. 
(S/N) 

1 

1A 

Nozzle  Se^nent  B,  Top 

PP81  (851  OB- 5) 

5 

4  (BA64) 

2 

IB 

Nozzle  Segment  B,  Center 

78HB  (8510B-5) 

5 

5  (BA68) 

3 

1C 

Nozzle  Segment  B,  Bottom 

79HB  (8510B-5) 

5 

6  (BA82) 

Roan  Air 


Vertical  Stilling  Chamber 


Horizontal  Stilling  Chamber 


Test  Section,  Top 


Test  Section,  Bottom 


Test  Section,  Row  1,  Col  1 


Test  Section,  Row  3,  Col  1 


Barometer 


89TA  (8510B-15) 


23LG  (8510B-100) 


74BF  (8S06B-5) 


PP67  (8510B-5) 


Test  Section,  Row  3,  Col  3 


Test  Section,  Row  1,  Col  4 


Test  Section,  Row  3,  Col  4 


Nozzle  Segment  A,  Top 


Nozzle  Segment  A,  Center 


Nozzle  Segment  A,  Bottom 


Nozzle  Segment  C,  Top 


Nozzle  Segment  C,  Center 


Nozzle  Segment  C,  Bottom 


97BF  (8506-5) 


79BF  (85068-5) 


M038  (8506A-5) 


Air  Thermocouple 


He  Thermocouple 


Air  Orifice  Meter,  Upstream 


He  Venturi  Meter,  Upstream 


Microstar  (Left) 


Microstar  (Right) 


WL44  (8530A-100) 


43YB  (8510B-500) 


10  (AE82) 


16  (BA57) 


11 

20 

Test  Section,  Row  1,  Col  2 

83BF  (8506B-5) 

5 

3  (AF03) 

12 

3A 

Test  Section,  Row  2,  Col  2 

HE99  (8506-5) 

5 

11  (AM03) 

13 

3B 

Test  Section,  Row  3,  Col  2 

75BF  (8506B-5) 

5 

15  (BA56) 

14 

3C 

Test  Section,  Row  1,  Col  3 

92BF  (8506B-5) 

5 

2  (AE77) 

15 

30 

Test  Section,  Row  2,  Col  3 

KL52  (8506-2) 

2 

12  (AE80) 

14  (AF84) 


7  (BA85) 


17  (BA55) 


25 

40 

He  Venturi  Meter,  Ownstrm 

Dynisco  Oifferential 

0-25psid 

20  (2310) 

26 

He  Venturi  Meter,  Upstream 

TCPT69-25 

6  V 

Vi  shay 

27 

5A 

Air  Orifice  Meter,  Upstream 

Statham  Oifferential 

t  1  psid 

0  (8A63) 

28 

Air  Orifice  Meter,  Ownstrm 

1472  P96-10-350 

10  V  max 

32 


5C 


100  psia 
500 


1  (BA67) 
18  (BA96) 


mass  flow  and  helium  mass  flow  rate  calculations  through  a 
square  edge  orifice  meter  and  venturi  meter.  The  pressure 
of  the  air  flow  at  the  upstream  tap  of  the  orifice  plate  was 
measured  using  an  ENDEVCO  Model  8530A-100.  The  differential 
pressure  across  the  orifice  plate  was  monitored  by  a  Statham 
Model  1472  P96-1D-350  differential  pressure  transducer  with 
a  ±1  psid  range.  Upstream  of  the  venturi  meter  throat  the 
helixim  pressure  was  measured  by  an  ENDEVCO  Model  8510 B- 500 
transducer.  The  differential  pressure  across  the  venturi 
was  measured  using  a  Dynisco  differential  pressure 
transducer,  Model  TCPT69-25,  with  a  0  to  25  psid  range. 

Signal  Conditioning  Units/Power  Supplies.  Five  ENDEVCO 
Model  4225  ac-operated  power  supplies  were  used  to  provide 
±18  volts  to  the  18  ENDEVCO  Model  4423  Signal  Conditioning 
Units  (SCUs) .  These  SCUs  provide  10  volts  excitation  to  the 
17  ENDEVCO  pressure  transducers  and  the  Statham  differential 
pressure  transducer  (position  27-28) .  The  SCUs  also  amplify 
output  voltage  from  each  transducer.  A  gain  of  20  was 
selected  for  these  transducers.  The  Dynisco  differential 
transducer  (position  25-26)  required  6  volts  excitation 
voltage  and  high  gain  in  order  to  resolve  the  low  pressure 
differential.  A  Vishay  Measurements  Group  Signal 
Conditioner  Amplifier,  Model  2310,  was  utilized  with  a  gain 
setting  of  500.  Table  3.3  lists  the  serial  numbers  of  each 
SCU  and  which  transducer  it  was  assigned  to. 
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Col  1 


Col  2 


Col  3 


Col  4 


Nozzle  Segment  A 


Nozzle  Segment  B 


Nozzle  Segment  C 


Figure  3.8  Test  Section  Transducer  Tap  Locations 


Nicolet  System  500.  The  Nicolet  System  500  is  a 
complete  Data  Acquisition  System.  A  pedestal  style  Data 
Acquisition  Unit  (part  no.  886-501)  houses  one  540  CPU  board 
and  five  digitizer  boards.  Each  of  these  Nic  514  digitizer 
boards  has  a  four  channel  capability  making  a  total  of  20 
channels  available  for  simultaneous  sampling  of  analog 
signals.  Once  digitized,  signals  are  transferred  through  a 
General  Purpose  Interface  Bus  (GPIB)  to  the  system 
controller.  The  controller  is  an  IBM  compatible  386 
computer  with  2  Megabytes  of  RAM,  a  150  Megabyte  hard  drive. 
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one  1.2  Megabyte  5.25  inch  floppy  disk  drive,  one  3.5  inch 
disk  drive,  and  a  VGA  video  display  terminal.  Nicolet 
system  software,  version  5.10,  runs  in  a  Microsoft 
Windows/386  format  to  provide  control  of  all  aspects  of  the 
acquisition  and  display  of  data.  For  air-only  runs 
triggering  of  the  Nicolet  System  500  occurred  when  voltage 
output  of  the  transducer  placed  in  the  horizontal  stilling 
chamber  (position  6)  began  to  rise  from  its  initial  steady- 
state  condition  indicating  opening  of  the  3-inch  ball  valve 
in  the  air  supply.  The  triggering  criteria  for  air-helium 
runs  was  the  voltage  level  output  from  the  helium  venturi 
meter  differential  pressure  transducer  (position  25-26) .  A 
step  increase  in  voltage  from  this  transducer  indicated 
onset  of  helium  flow  and  started  data  collection. 

Two  computer  programs,  written  in  GW-Basic,  control  the 
pre-run  inputs  and  the  post-run  data  reduction.  Appendix  B 
contains  the  program  listings. 

Schlieren  Optical  System 

Schlieren  photography  and  filming  of  the  test  section 
flow  dynamics  were  accomplished  using  the  setup  shown  in 
Figure  3.9.  A  Gordin  spark  lamp.  Model  5401,  provided 
lighting  for  still  photographs  using  Polaroid  3000  ASA  film. 
A  Gordin  Power  Supply,  Model  5205,  powered  the  spark  lamp. 
The  two  mirrors  were  7.5  inch  diameter  concave  with  40  inch 
focal  lengths.  A  second  set  of  10  inch  diameter,  concave. 
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30  inch  focal  length  mirrors  was  also  used.  A  steady  300 
watt  arc-lamp  point  source  illuminated  the  test  section  for 
filming.  A  Redlake  16mm  intermittent  high  speed  motion 
picture  camera.  Model  51-0003,  Serial  Number  647,  was  placed 
just  behind  the  knife-edge  and  the  image  focused  on  the  film 
plane.  Film  speeds  ranged  from  100  to  500  frames  per 
second. 


Figure  3.9  Schlieren  Optical  System 
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IV.  Experimental  Procedures 


The  experimental  procedures  followed  during  this  research 
included  calibration  of  pressure  transducers,  schlieren 
photography  techniques,  and  data  acquisition  and  reduction 
steps. 

Pressure  Transducer  Calibration 

AMETEK  PneuE. 'tic  Dead  Weight  Testers,  Model  PK-2  (low 
pressi  re,  300  inches  of  water)  and  Model  HK-500  (high 
pressure,  500  psig) ,  were  used  to  calibrate  all  transducers. 
The  manufacturer  of  these  pneumatic  dead  weight  testers 
cautions  that  the  use  of  any  gas  other  than  clean  compressed 
air  or  nitrogen  could  cause  discrepancies  in  the  calibration 
data  (10) .  The  pressure  transducers  were  calibrated  with 
their  respective  SCUs  in  line.  Since  the  SCU  does  not  have 
a  unity  gain  setting,  it  will  always  multiply  the  pressure 
transducer  output  voltage  by  some  factor  depending  on  which 
gain  setting  has  been  selected.  The  ENDEVCO  SCUs  have  four 
gain  settings  (5,  10,  20,  50)  and  a  bridge  balance 
adjustment  screw.  The  Vishay  2310  Signal  Conditioner 
Amplifier  has  a  variable  gain  adjustment  and  was  set  at  500. 
The  gain  factors  were  not  exactly  as  indicated  on  the  front 
panel;  therefore,  the  pressure  transducers  were  calibrated 
with  the  SCU  in  line  so  that  the  exact  gain  factor  is 
included  in  the  slope  of  the  calibration  curve  (4:4-1).  The 
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calibration  curves  of  the  pressure  transducers  used  in  this 
research  are  presented  in  Appendix  A. 

Schlieren  Photographic  Technicmea 

The  setup  of  the  schlieren  equipment  shown  in  Figure  3.9  was 
accomplished  using  the  guidelines  of  AGARDograph  No.  23 
(11:13-15).  Obtaining  quality  photographs  involved  much 
trial  and  error  due  to  the  extremely  low  density  and  small 
optical  depth  in  the  test  section.  For  high  speed  filming, 
it  was  discovered  that  placing  the  camera  directly  'renind 
the  knife-edge's  focal  point  provided  the  best  focus  and 
lighting  capabilities.  An  attempt  to  film  directly  off  the 
image  projected  on  the  still  camera's  ground  glass  captured 
some  flow  dynamics;  however,  the  lighting  was  severely 
diminished.  Various  film  speeds  were  examined  to  determine 
if  an  optimum  existed.  Results  were  very  similar  for  film 
speeds  in  the  range  of  100  to  500  frames  per  second. 

Data  Acquisition  and  Reduction 

Figure  4.1  shows  the  steps  used  in  acquiring  and  reducing 
data.  Details  concerning  the  use  of  Nicolet's  System  500 
are  contained  in  the  user's  operation  manual  (12) .  The 
Nicolet  software  runs  in  a  Windows  3.0  environment.  For 
further  information  on  Windows  procedures  refer  to  the 
manual  (13) . 


4 


2 


Figure  4.1  Data  Acquisition  Procedures 
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V.  Results  and  Discussion 


Results  of  the  research  accomplished  include  evaluation 
of  the  simulated  laser  cavity  test  equipment,  ex2uiiination  of 
the  schlieren  still  and  high  speed  photography,  and  analysis 
of  pressure  data  obtained  during  runs  made  under  various 
flow  conditions. 

Test  Ecmipment  Performance 

Monitoring  of  pressure  data  over  a  one  hour  period 
revealed  the  stable  characteristics  of  the  pressure-vacuum 
flow  system.  The  100  psig  air  pressure  supply  exhibited  a 
cyclic  fluctuation  of  approximately  5  minutes  period  and  2.5 
psi  amplitude.  This  cycle  appeared  to  correspond  to  the 
on/off  cycling  of  the  air  dryers  connected  to  the  air  supply 
line.  The  period  was  not  considered  detrimental  as  test 
runs  were  only  of  10  seconds  duration  and  the  air  supply 
pressure  during  operation  was  observed  to  remain  steady. 

The  vacuum  pumps  were  capable  of  achieving  a  minimum 
pressure  inside  the  cavity  of  approximately  6  mm  Hg  (0.116 
psia)  in  roughly  25  minutes  of  evacuation.  Cavity  pressure 
was  monitored  via  a  calibrated  vacuum  gauge  prior  to  each 
run  to  determine  a  steady  condition  and  maximum  vacuum. 

The  reference  vacuum  applied  to  each  psig  pressure 
transducer  was  monitored  by  a  calibrated  portable  vacuum 
standard.  The  minimum  reference  pressure  achieved  was 
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stable  at  1.0  mm  Hg  (0.0193  psia)  in  approximately  30 
minutes  of  pumping.  The  selection  of  transducers  with 
reference  ports  that  did  not  leak  was  essential  to  maintain 
the  lowest  possible  reference  vacuum.  ENDEVCO  transducer 
reference  ports  were  tested  for  leakage  using  a  hand  pumped 
vacuum  gauge  prior  to  selection  for  use. 

Early  testing  of  the  helitim  supply  line  disclosed  a 
tendency  for  fluctuations  from  the  air  supply  to  travel 
upstream  to  the  helium  venturi  flow  meter  which  caused 
erroneous  mass  flow  rate  calculations.  To  overcome  this,  a 
helium  flow  control  nozzle  was  designed  to  choke  the  stream 
prior  to  entering  the  horizontal  stilling  chamber.  Mass 
flow  of  the  helium  was  controlled  by  the  second  dome  valve 
pressure  on  the  helium  supply  line.  This  pressure  (P^^) , 
set  at  the  helium  farm  regulator  board,  was  varied  from  zero 
for  air-only  test  runs  to  100  to  350  psig  for  runs  with 
helium  injected  into  the  mixer.  Data  runs  made  with  the 
supply  pressure  set  below  approximately  175  psig  failed  to 
prevent  the  fluctuations  from  travelling  upstream  as  this 
was  the  pressure  where  choking  of  the  flow  control  nozzle 
occurred.  Data  for  these  runs  was  examined;  but,  was  not 
considered  reliable  enough  for  use. 
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Schlieren  Photography  Observatlona 

Flow  dynamics  inside  the  simulated  laser  cavity  were 
viewed  at  single  time  points  using  still  schlieren 
photography.  High  speed  schlieren  filming  also  captured  the 
flowfield  for  the  duration  of  ten  second  runs.  The  quality 
of  these  photographic  endeavors  was  marginal.  It  is 
believed  that  the  very  low  pressures,  and  consequently  the 
low  density  gradients,  in  the  test  cavity  as  well  as  the 
poor  cjuality  plexiglas  walls  contributed  to  the  inferior 
pictures.  Dark  zones  were  abundant  in  many  of  the 
photographs  and  obtaining  a  uniform  field  was  not  possible. 
Despite  this,  much  information  was  gleaned  from  the  many 
photographs  and  film  runs  made.  Because  the  reproducibility 
of  the  photographs  was  limited,  presentation  of  these 
observations  will  be  made  schematically. 

Figure  5.1  illustrates  the  steady  state  flowfield  in 
the  cavity  which  occurred  for  all  visualized  runs  during  at 
least  the  first  2  seconds  of  operation.  Jet  boundaries 
develop  at  the  exit  of  the  nozzles  forming  triangular  base 
zones  in  which  recirculation  of  the  flow  can  take  place.  At 
the  apex  of  each  triangle  a  base  wake  region  appears 
followed  by  oblique  shocks.  Varying  the  value  changed 
the  conditions  beyond  the  2  second  range.  These  effects 
will  be  discussed  in  the  next  section  concerning  pressure 
data  results. 
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Figure  5.1  Sche*^®*'*®  Cavity  Flow  at  Approximately  2  Seconds 


The  area  observed  in  filming  included  a  circular 
portion  of  the  cavity  from  the  two  whole  nozzle  exits 
downstream  to  approximately  half  the  test  section  length. 
Figure  5.2  depicts  the  flow  dyneu&ics  from  a  typical  10 
second  high  speed  film  run.  The  flow  rapidly  established 
the  pattern  of  Figure  5.1  within  the  first  second  of  oper¬ 
ation  and  maintained  this  wave  formation  for  approximately  2 
seconds.  As  the  cavity  pressure  increased  during  a  test  run 
density  gradients  became  larger  improving  the  visibility  of 
the  jet  boundaries  and  oblique  shocks.  These  jets  and 
shocks  were  also  observed  in  still  schlieren  photographs. 
Following  the  2  seconds  of  apparent  steady  flow,  transient 
wave  pulses  were  observed  exiting  the  nozzles  causing 
disturbance  of  the  previously  stable  oblique  shocks.  From 
approximately  3  seconds  to  6  seconds  these  transients 
continued  while  the  uppermost  oblique  shock  began  rotating 
into  the  cavity.  The  shock  could  be  seen  beginning  to 
nearly  impinge  on  the  transducer  at  position  14  at  approx¬ 
imately  5.0  seconds  after  flow  initiation.  It  is  possible 
that  boundary  layer  growth  along  the  cavity  ceiling  builds 
up  as  the  pressure  in  the  cavity  begins  to  rise  forcing  the 
oblique  shock  downward  at  a  sharper  angle.  As  the  vacuum 
system  begins  to  fill,  the  back  pressure  rises  and  the  base 
wake  begins  to  open  (t»7.0  seconds).  The  wake  becomes  fully 
open  at  approximately  10  seconds. 
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Pressure  Data 


Voltage  readings  from  piezoresistive  pressure 
transducers  placed  in  the  test  apparatus  were  converted  to 
pressure  values  by  the  computer  progreun  POSTRUN. BAS  listed 
in  Appendix  B.  This  program  also  calculated  mass  flow  rates 
of  the  air  and  helium  supplies,  nozzle  exit  Mach  numbers  and 
mixture  specific  heat  ratios  during  a  run. 

Table  5.1  lists  the  run  ntimbers  and  operating 
conditions  for  the  most  significant  runs  made.  These  runs 
were  selected  as  most  illustrative  of  effects  of  changes  in 
specific  heat  ratio  (y) ,  and  effects  of  the  presence  or 
absence  of  the  air/helium  mixer. 


Table  5.1  Run  Numbers  and  Operating  Conditions 


Run  Nvimber 

p„,  (psig) 

Mixer 

020 

0  Air  Only 

Yes 

023 

175 

Yes 

024 

200 

Yes 

025 

225 

Yes 

026 

250 

Yes 

027 

275 

Yes 

028 

300 

Yes 

029 

350 

Yes 

034 

0  Air  Only 

No 

035 

225 

No 
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To  control  the  mixture  y  value  in  the  cavity  the  mass 
flow  rate  of  helium  was  varied  by  increasing  the  helium 
injection  pressure,  Figure  5.3  shows  the  effect  of 

increasing  on  the  mixture  y  value.  Trends  are 
illustrated  using  a  seventh  degree  polynomial  best  fit  curve 
through  all  data  points  for  a  given  P^^.  The  design  y  value 
of  1.51  is  approximately  achieved  by  using  a  P^^  pressure  of 
200  psig.  Increasing  the  helium  injection  pressure 
increases  the  helium  mass  flow  which  consequently  raises  the 
Y  value. 

Figure  5.4  shows  calculated  (experimental)  nozzle 
exit  Mach  numbers  for  various  helium  injection  pressures. 
Trends  again  are  depicted  using  a  seventh  degree  polynomial 
best  fit  curve  through  all  data  points;  however,  to  retain 
clarity  only  one  third  of  the  data  points  are  displayed.  The 
flow  readily  established  velocities  in  the  hypersonic  range. 
Higher  nozzle  exit  Mach  numbers  are  achieved  for  increased 
Pjj^  values.  This  is  attributed  to  the  increased  gas  y 
values. 

When  results  of  experimental  and  theoretical  exit  Mach 
numbers  are  compared  (equations  7  and  9,  chapter  2)  a 
Reynolds  number  effect  is  discovered.  Figure  5.5  plots 
experimental  and  theoretical  nozzle  exit  Mach  numbers  for  a 
range  of  y  from  1.4  (air  only)  to  1.58  (Pg^=350  psig).  The 
boundary  layer  thickness,  ft,  is  proportional  to  the  inverse 
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Figure  5.3  Specific  Heat  Ratios  vs  Time  for  Various  Values 
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Nozzle  Exit  Mach  Number 


Seconds 


Figure  5.4  Nozzle  Exit  Mach  Numbers  for  Various  Values 
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Figure  5.5  Experimental  and  Theoretical  Nozzle  Exit  Mach 
Number  Comparison 


of  the  Reynolds  number  usually  by  a  power  law.  Reynolds 
numbers  at  the  nozzle  throat  were  approximately  1.4  x  10^ 
for  flow  with  air  only  (Re=Vh/v)  which  relates  to  thin 
boundary  layers  and  the  isentropic  assumption  is  valid  (14) . 
However,  as  helium  is  introduced  into  the  flow  y  increases 
as  well  as  the  acoustic  velocity  and  kinematic  viscosity  of 
the  flow  medium.  Also,  the  density  and  Reynolds  numbers 
decrease  which  relates  to  an  increase  in  6  along  the  nozzle 
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wall.  This  consequently  decreases  the  effective  area  ratio 
(e)  of  the  nozzle  and  the  exit  Mach  number  diverges  from  the 
isentropic  case  as  seen  in  Figure  5.5. 

Figure  5.6  plots  the  horizontal  stilling  chamber 
pressure  during  the  first  3  seconds  of  operation  for  various 
values.  Pressures  remain  at  their  3  second  value  for 
the  remainder  of  the  10  second  run.  Three  runs  in  which 
helium  was  not  injected  into  the  system  (P„,=0)  demonstrate 
the  repeatability  of  operating  the  test  apparatus.  Repeat 
runs  were  made  for  each  condition  examined  with  similar 
results;  however,  these  were  .lot  plotted  to  retain  clarity 
in  the  graphs.  The  air-helium  runs  indicated  a  slightly 
increased  chamber  pressure  for  increasing  P^^  settings. 

As  helium  is  introduced  into  the  mixer  there  is  an  increased 
lag  before  the  initial  pressure  rise  in  the  horizontal 
stilling  chamber.  This  is  probably  due  to  the  different 
methods  used  in  triggering  the  data  acquisition  sequence  for 
air-only  or  air-helium  runs  (refer  to  page  3-17). 

To  compare  cavity  pressure  conditions  under  various  P^^ 
levels  graphs  of  pressure  versus  time  for  several  of  the 
runs  listed  in  Table  5.1  are  presented  in  Figures  5.7 
through  5.15.  Each  graph  contains  recorded  pressure  values 
in  the  lasing  cavity  over  a  ten  second  run  viewed  as  a 
population.  Stable  pressure  conditions  in  the  cavity  are 
characterized  as  low  level  (below  approximately  0.25  psia) 


5 


12 


Horizontal  Stilling  Chamber  Pressure  (psia) 


110 
100 
90 
80 
70 
60 
50 
40 
30 
20 
10 
0 

0.0  0.5  1.0  1.5  2.0  2.5  3.0 

Seconds 

Figure  5.5  Horizontal  Stilling  Chamber  Pressure  vs  Time  for 
Various  Settings 
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maintained  for  as  long  a  possible. 

Figure  5.7  represents  data  obtained  during  run  029 
which  used  a  helixim  inlet  pressure  of  350  psig  achieving  an 
approximate  y  of  1.58.  This  condition  was  considered  the 
upper  operating  limit  of  the  test  apparatus  as  unstable 
pressures  prevailed  and  no  stable  performance  time  was 
observed.  This  condition  was  compared  to  other  operating 
conditions  but  will  not  be  presented  in  any  further  data 
correlations  as  no  beneficial  information  was  found. 
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Figure  5.7  Cavity  Pressures  vs  Time  for  Run  029 
(P„.=350  psig,  y«1.58) 
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To  examine  the  opposite  end  of  the  helitim  injection 
pressure  range,  a  run  was  done  with  air  only.  Figure  5.8 
plots  pressures  in  the  cavity  during  run  020,  a  typical  ten 
second  air-only  rvin.  Approximately  5  seconds  of  stable 
pressure  performance  occurs  before  the  pressure  fluctuations 
seen  in  Figure  5.8  appear.  A  pressure  peak  of  approximately 
1.0  psia  at  0.1  second  could  indicate  a  shock  passing 
through  the  system  at  start-up.  Another  around  4.0  seconds 
could  be  a  precursor  to  the  unstable  behavior  at  5.0 
seconds.  This  would  seem  to  indicate  that  the  process  of 
injecting  helium  into  the  mixer  was  at  least  partly 


Seconds 


Figure  5.8  Cavity  Pressures  vs  Time  for  Rvin  020 
(P^^=»0,  Air  Only,  y»1.4) 
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responsible  for  the  onset  of  the  unfavoreUsle  pressure 
behavior. 

To  discover  another  possible  factor  which  may  have 
caused  the  pressure  fluctuations  to  occur,  the  mixer  was 
removed  and  a  ten  second  run  was  accomplished  using  the  same 
air-only  condition.  Run  034  cavity  pressure  measurements 
are  displayed  in  Figure  5.9.  Stable  performance  occurred 
for  the  duration  of  the  run  which  means  that  the  mixer  was 
also  contributing  to  the  onset  of  fluctuations  as  seen  in 
Figure  5.8. 
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Figure  5.9  Cavity  Pressures  vs  Time  for  Run  034 

Air  Only,  yfsl.A,  No  Mixer) 
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So  stable  run  tine  is  affected  by  the  helium  injection 
pressure  and  by  the  use  of  the  mixer.  Stable  performance 
should  be  seen  for  runs  made  with  settings  between  the 
extremes  of  0  and  350  psia.  Figures  5.10  through  5.13 
verify  this  trend  for  runs  made  with  the  mixer  installed. 

Figure  5.10  shows  pressure  data  in  the  cavity  during 
run  023  which  used  a  value  of  175  psig  corresponding  to 
an  approximate  y  value  of  1.48  (see  Figure  5.3).  There  was 
approximately  2 . 5  seconds  of  stable  performance  for  this 
run. 
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Figure  5.10 
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Cavity  Pressures  vs  Time  for  Run  023 
(P^*175  psig,  y»1.48) 
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Run  024  pressure  data  (Figure  5.11)  for  a  setting 
of  200  psig  reveals  approximately  the  same  stable  time  as 
run  023,  2.5  seconds.  The  25  psig  delta  from  175  to  200 
psig  did  not  influence  the  stability  time.  It  should  be 
noted  that  this  case  provided  the  y  value  most  near  the 
nozzle  design  y  of  1.51. 
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Figure  5.11  Cavity  Pressures  vs  Time  for  Run  024 
(P^,=200  psig, 


Increasing  the  helium  injection  pressure  to  225  psig 
(YsI.54)  decreased  the  stable  rxin  time  to  approximately  2.25 
seconds  as  seen  in  Figure  5.12. 
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Figure  5.12  Cavity  Pressures  vs  Tine  for  Run  025 
(P^^=225  psig,  Y*i‘54) 

Figure  5. 13  displays  the  decrease  in  stable  run  time  to 
1.8  seconds  as  P^^  was  increased  to  275  psig  in  run  027 
(Y«1.56)  . 

A  final  configuration  was  examined  in  which  the  mixer 
was  removed  and  helium  was  injected  directly  into  the  air 
stream.  Mixing  of  the  air  and  helium  was  then  by  turbulent 
and  molecular  diffusion.  Figure  5.14  shows  the  pressure 
data  for  this  case,  run  035,  which  used  a  P^^  of  225  psig. 
When  compared  to  run  025  (Figure  5.12)  the  influence  of  the 
mixer  on  the  onset  of  the  pressure  fluctuations  becomes 
evident.  An  increase  in  stable  run  time  to  approximately 
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Figure  5.13  Cavity  Pressures  vs  Time  for  Run  027 
(P^^=275  psig,  Y»1.56) 
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Cavity  Pressures  vs  Time  for  Run  035 
(P^^*225  psig,  Nc  Mixer) 
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6.0  seconds  occurred. 

This  would  indicate  that  a  different  form  of  mixing  the 
air  and  helium  supplies  is  warranted  if  test  times  with 
longer  stability  are  desired.  However,  it  should  be 
understood  that  a  stable  run  time  of  approximately  2  seconds 
is  long  compared  to  the  short  duration  of  most  chemical 
lasers.  Reference  15  gives  pulse  durations  in  the 
microsecond  range. 

A  closer  examination  of  the  y  values  plotted  in  Figure 
5.3  may  reveal  a  relation  between  the  oscillating  trends  of 
those  ys  and  the  pressure  fluctuations  in  the  cavity. 

Figure  5.15  graphs  y  during  a  3  second  duration  for  3 
cases,  Pj,^=200  psig,  P„^=225  psig,  and  P^g*225  psig  sans 
mixer.  Note  that  the  case  in  which  the  mixer  was  removed 
does  not  provide  a  stable  y  as  it  begins  lower  than  its 
counterpart  with  the  mixer  installed  and  slowly  increases  to 
approximately  1.53.  Also  notice  the  oscillating  pattern  of 
the  trends  represented  using  a  seventh  degree  polynomial 
best  fit  curve.  Since  the  values  for  y  are  calculated  from 
experimental  data,  in  particular  the  mass  flow  rates  of  air 
and  helium  for  a  given  run,  an  exeunination  of  those 
variables  is  warranted. 

Mass  flow  rates  for  rvms  025  and  035  which  used  a  P„^ 
setting  of  225  psig  with  and  without  the  mixer  installed  are 
plotted  in  Figure  5.16.  The  mass  flow  rate  for  air  in  both 
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Figure  5.15  Specific  Heat  Ratio  Comparison 
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Figure  5.16  Air  and  Helium  Mass  Flow  Rates 
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cases  remains  fairly  stable  with  a  slightly  negative  slope. 
The  helivim  mass  flow  rates  however  exhibit  trends  similar  to 
those  of  Y  in  Figure  5.15.  Note  the  same  increasing  trend 
in  the  helium  mass  flow  for  the  non-mixer  case  as  seen  in 
the  Y  value. 

The  mass  flow  calculation  for  helium  is  a  function  of 
the  helium  supply  pressure  upstream  of  the  venturi  meter  and 
the  differential  pressure  across  the  venturi  throat. 
Examination  of  the  helium  supply  pressure  (position  32) 
revealed  stable  pressure  at  approximately  500  psia.  When 
the  differential  pressure  across  the  venturi  (positions  25- 
26)  was  plotted  (Figure  5.17)  the  trends  were  nearly 
identical  to  that  3f  the  y  values  (Figure  5.15).  Further 
investigation  of  the  helium  supply  pressure  and  differential 
pressure  is  needed  to  determine  whether  a  direct  relation 
exists  between  the  oscillations  seen  in  y  trends  and  the 
pressure  fluctuations  seen  in  the  cavity. 

Focussing  attention  on  the  initial  pressure  responses 
in  the  cavity  during  the  first  two  seconds  of  operation,  the 
following  graphs  of  pressure  versus  time  at  each  transducer 
location  were  generated.  Figures  5.18  through  5.30  display 
the  pressure  history  at  each  position  for  air-helium  runs 
023,  024,  025,  and  027  made  with  the  mixer  installed  to 
determine  possible  effects  of  varying  y*  Refer  to  Figure 
3.8  for  position  number  locations. 
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Figure  5.18  displays  static  pressure  histories  for 
position  1.  Behavior  for  each  run  is  similar  and  no  effect 
of  increasing  y  can  be  discerned.  An  initial  pressure  drop 
to  approximately  0.04  psia  indicates  the  starting  of  the 
system  followed  by  a  gradual  increase  in  pressure  to  a 
stable  level  at  approximately  0.1  psia.  This  trend  of 
gradually  increasing  pressure  is  a  reflection  of  the 
upstream  horizontal  stilling  chamber  pressure  which  is  in 
transition  to  its  maximum  as  previously  seen  in  Figure  5.6. 

The  base  pressure  at  position  2  (Figure  5.19)  also 
drops  initially  and  increases  gradually.  The  pressure  in 
this  base  region  fluctuates  more  than  at  position  1  which 
could  be  attributed  to  the  recirculation  characteristic  of 
the  flow  at  the  base. 

Position  3  pressures  shown  in  Figure  5.20  resemble 
those  recorded  at  position  1  with  no  apparent  y  effect 
observed. 

Nozzle  exit  pressures  at  positions  17  and  18  are 
plotted  in  Figures  5.21  and  5.22  respectively.  Behavior 
reflects  that  seen  from  the  base  transducers  (positions  1 
through  3)  and  y  seems  to  have  no  effect  on  the  amount  of 
initial  pressure  drop. 

Moving  downstream,  into  the  cavity,  2  second  pressure 
histories  for  column  3  of  the  test  cavity  are  presented  in 
Figures  5.23  through  5.25. 
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Figure  5,18  Position  1  Pressure  History  for  Various  P^^  Values 


Figure  5.19  Position  2  Pressure  History  for  Various  P^  Values 
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PSIA  :  PSIA 


>.20  Position  3  Pressure  History  for  Various  P^^  Values 


Figure  5.21  Position  17  Pressure  History  for  Various  P^  Values 


Seconds 

Figure  5,22  Position  18  Pressure  History  for  Various  P^^  Values 


Position  14  (Figure  5.23)  pressures  are  similar  to 
those  observed  upstream.  This  is  contrary  to  what  is 
expected  from  theory.  If  the  flow  is  expanding  from  the 
nozzle  exit  into  the  cavity,  pressure  should  decrease 
downstream  relating  to  higher  velocity.  This  phenomenon 
will  be  addressed  further  in  a  later  discussion  about  Mach 
number.  The  onset  of  the  pressure  fluctuations  is  seen  at 
approximately  1.6  seconds  for  run  027  when  P^^  is  set  at  275 
psig. 

Figure  5.24  displays  the  pressures  for  position  15. 
After  the  initial  drop,  pressures  rapidly  increase  as  the 
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Figure  5.25  Position  16  Pressure  History  for  Various  P^^  Values 

base  wake  is  established  and  the  oblique  compression  shocks 
form.  The  onset  of  the  unstable  pressure  fluctuations  is 
again  seen  at  approximately  1.6  seconds  for  the  275  psig 
case. 

An  interesting  event  is  observed  in  the  lower  row  of 
the  cavity  as  seen  in  the  pressure  trends  of  position  16 
(Figure  5.25).  The  runs  with  P^^  set  at  200,  and  225  psig 
have  a  pressure  spike  at  0.3  seconds.  Since  the  magnitude 
of  the  spike  does  not  correspond  to  an  increase  in  y  the 
effect  is  probably  not  attributable  to  that  variable. 

Instead,  the  presence  of  a  weak  shock  passing  through  the 
cavity  during  start-up  is  a  more  probable  cause.  This 
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start-up  shock  may  not  be  captured  during  data  acquisition 
due  to  its  short  duration  and  therefore  does  not  appear  in 
other  run's  pressure  results. 

Figures  5.26  through  5.28  show  pressures  for  various 
settings  downstream  along  column  2  in  the  cavity. 

Increased  pressure  trends  in  Figure  5.26  would  indicate 
that  the  oblique  shocks  shown  in  Figure  5.1  may  actually  be 
crossing  near  position  11. 

Position  12  pressures  depicted  in  Figure  5.27  reflect 
reasonably  stable  performance  with  a  dip  in  pressure  of 
approximately  0.05  psia  followed  by  a  gradual  increase. 
Boundary  layer  buildup  along  the  cavity  walls  could  be 
acting  to  diffuse  the  flow;  lowering  the  velocity  and 
raising  the  pressure. 

Pressure  data  at  position  13  in  row  3  of  the  cavity 
(Figure  5.28)  would  also  indicate  oblique  compression  shocks 
crossing  at  that  point. 

Position  9  and  10  pressures  are  shown  in  Figures  5.29 
and  5.30  respectively.  The  onset  of  unstable  pressure 
behavior  at  both  of  these  positions  is  again  observed  for 
the  275  psig  case  at  approximately  1.6  seconds.  Position  10 
pressures  shown  in  Figure  5.30  show  the  initial  pressure 
peak  previously  seen  which  could  be  caused  by  a  shock 
passing  through  the  cavity. 
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Figure  5.28  Position  13  Pressure  History  for  Various  P  Values 
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Figure  5.29  Position  9  Pressure  History  for  Various  P  Values 
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Figure  5,30  Position  10  Pressure  History  for  Various  P^^  Values 

To  further  examine  the  expected  expansion  of  the  flow 
into  the  cavity,  Mach  numbers  at  position  14  and  16  were 
calculated  assuming  isentropic  flow  to  those  locations. 

These  Mach  numbers  are  compared  to  Mach  numbers  calculated 
at  the  nozzle  exits  in  Figures  5.31  and  5.32  for  the  cases 
of  Pj,^  at  225  psig  with  and  without  the  mixer  installed. 

With  the  mixer  present  (Figure  5.31),  unsymmetrical 
behavior  in  the  cavity  was  observed  as  the  row  1  Mach  number 
at  position  14  was  greater  than  the  Mach  number  at  position 
17  indicating  expansion  as  the  flow  moved  downstream.  In 
the  lower  portion  of  the  cavity,  row  3  Mach  numbers  remained 
fairly  constant. 
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Mach  Number  Mach  Number 
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Figure  5.32  Cavity  Mach  Number  Comparison  sans  Mixer 
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Mach  numbers  at  position  14  and  16  with  the  mixer 
removed  reflected  their  upstream  levels  as  no  distinct 
indication  of  expansion  is  seen  in  Figure  5.32. 

Flow  dynamics  in  the  test  cavity  are  very  complex. 
Boundary  layer  growth  along  the  upper  and  lower  walls  of  the 
cavity  as  well  as  along  the  vertical  sidewalls  can 
contribute  to  increasing  3  dimensional  effects  as  the  flow 
moves  downstream.  This  can  cause  a  diffusing  effect  which 
cancels  out  the  expected  expansion  into  the  cavity. 
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VI .  Conclusions 


An  experimental  investigation  of  the  effect  of  varying 
the  specific  heat  ratio  for  the  flow  through  a  simulated 
chemical  laser  cavity  was  accomplished.  A  survey  of  static 
pressures  in  the  cavity,  still  schlieren  photographs,  and 
high  speed  schlieren  film  of  the  flowfield  determined  the 
behavior  of  the  flow.  The  following  conclusions  were 
reached: 

(a)  The  equipment  designed  to  supply  helium  to  the  system 
functioned  to  provide  control  of  the  mixture  specific 
heat  ratio  of  the  air-helium  flow.  However,  it  was  not 
determined  if  the  injection  of  helium  into  the  system 
was  the  direct  cause  of  unstable  pressure  behavior. 

(b)  For  specific  heat  ratios  ranging  from  approximately 
1.40  to  1.56  the  experimental  apparatus  provided 
similar  cavity  flow  conditions.  The  length  of  stable 
run  time  (before  pressure  instabilities  began) 
decreased  as  specific  heat  ratio  was  increased.  The 
upper  limit  on  specific  heat  ratio  for  this  system  was 
approximately  1.58  corresponding  to  a  helium  injection 
pressure  of  350  psig  which  gave  no  favorable  pressure 
behavior.  At  the  nozzle  design  specific  heat  ratio  of 
1.51  the  stable  run  time  was  approximately  2.5  seconds. 
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(c)  Increased  specific  heat  ratio  related  to  decrease 
Reynolds  number  and  consequently  increased  boundary 
layer  thickness.  This  caused  a  decrease  in  effective 
area  ratio  of  the  nozzles,  thus  exit  Mach  numbers  were 
lower  than  the  isentropic  assumption  predicted  as 
specific  heat  ratio  was  raised  from  the  level  of 
approximately  1.4. 

(d)  Boundary  layer  effects  along  the  sidewalls  of  the  test 
section  tended  to  cancel  out  pressure  expansion  in  the 
cavity.  However,  unsymmetrical  pressure  behavior 
indicated  some  expansion  in  the  upper  half  of  the 
cavity  when  the  mixer  was  installed. 

(e)  Use  of  the  mixer  decreased  stable  run  time  by 
introducing  unstable  pressure  fluctuations  earlier. 
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VII.  peeQmTOftndationa 


Further  investigation  of  the  cavity  flow  dynamics  in 
the  simulated  laser  is  recommended.  To  continue  this 
research  the  following  suggestions  are  presented; 

(a)  Since  it  is  desired  to  have  the  cavity  operate  at  the 
nozzle  design  specific  heat  ratio  of  1.51,  a  longer 
stable  testing  time  would  be  beneficial.  To  achieve 
this,  a  new  method  of  mixing  the  secondary  helium 
supply  with  the  primary  air  should  be  explored.  A  ring 
manifold  mounted  inside  the  10  inch  diameter  air  supply 
pipe  with  small  helium  supply  perforations  around  its 
circumference  could  provide  sufficient  mass  flow 
without  adverse  pressure  effects.  Or  direct  injection 
of  the  secondary  helium  at  circumferential  stations 
through  the  pipe  itself  could  serve  as  a  delivery 
system  without  positioning  a  body  directly  in  the  flow 
of  primary  air.  Also,  the  addition  of  a  concentration 
probe  cou'.d  be  used  to  verify  the  mixture  specific  heat 
ratio. 

(b)  The  initial  seconds  of  flow  through  the  cavity  is 
important.  More  rapid  filling  of  the  horizontal 
stilling  chamber  upstream  of  the  supersonic  nozzles 
would  bring  the  system  to  operating  condition  more 
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readily.  By  decreasing  the  volume  of  the  system  prior 
to  the  test  cavity  the  existing  mass  flow  of  air  would 
facilitate  a  more  rapid  pressure  response.  Shortening 
of  the  horizontal  stilling  chamber  downstream  of  the 
mixer  location  would  benefit. 

(c)  Addition  of  hypersonic  wedges  to  more  closely  model  the 
actual  chemical  laser  nozzle  system.  Milling  of  the 
existing  nozzle  segments  to  attach  several  rows  of 
hypersonic  wedges  would  be  a  logical  step  in  furthering 
this  research. 

(d)  To  improve  flow  visualization,  the  plexiglas  walls  of 
the  test  cavity  should  be  replaced  with  optical  glass. 
The  superior  quality  should  greatly  improve  schlieren 
photography  and  high  speed  film  of  the  cavity  flow 
dynamics.  Also,  increasing  the  width  of  the  lasing 
cavity  from  its  present  value  would  provide  greater 
bending  of  the  light  as  a  result  of  the  density 
gradient  and  consequently  improve  the  flow  features 
captured  by  the  cameras. 
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The  calibration  curves  of  the  nineteen  pressure  transducers 
used  in  this  investigation  are  presented  in  the  following 
figures: 


Figure 

Number 


A.l  Calibration  Curve  for  Pressure  Transducer  Serial 
PP81 


Figure 

Number 


Figure 

Number 


A. 2  Calibration  Curve  for  Pressure  Transducer  Serial 
78HB 


A. 3  Calibration  Curve  for  Pressure  Transducer  Serial 
79HB 
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Figure  A. 10  Calibration  Curve  for  Pressure  Transducer  Serial 
Number  75BF 


Figure  A. 11  Calibration  Curve  for  Pressure  Transducer  Serial 
N\imber  92BF 
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Appendix  B:  Data  Acauialtlon  Software 


The  two  major  programs  which  control  pre-run 
acquisition  and  post-run  reduction  of  data  were  written  in 
Microsoft  GW-BASIC. 


PRERUN. BAS 


10  REM  Progran  PRERUN. BAS 
20  REM 

30  REM  Written  by:  Curt  0.  Botte 

40  REM  This  program  queries  the  user  for  pre-run  infonestion  necessary  for 
50  REM  post  data  acquisition  calculations.  This  prograa  should  be  run  and 
60  REM  information  provided  as  close  as  possible  to  actual  data  acquisition. 
70  REM  Each  time  this  program  is  run  a  COUNTER  file  is  accessed  and  the  rm 
80  REM  nunber  is  incremented.  An  output  file  containing  the  following  inputs 
90  REM  is  created  with  e  "PRE(run  reaper). DAT*  filerume  format.  This  output 
100  REM  file  will  be  accessed  during  the  POSTRUN  data  reduction. 

110  REN 

120  REN  1)  The  Heliue  Line  Tempereture  -  this  value  is  the  degrees 
130  REM  Fahrenheit  displayed  on  the  right  hand  Microstar  F 

140  REN  digital  display.  The  line  should  be  pressurized 

150  REN  with  heliuB  prior  to  recordirig. 

160  REM 

170  REM  2)  The  Air  Line  Teeperature  -  similarly,  this  value  is  read 
180  REM  off  the  left  hand  Microstar  F  digital  display.  The 

190  REN  air  line  is  always  pressurized  with  100  psis  air  at 

200  REN  the  thermocouple  location. 

210  REM 

220  REM  3)  Atmospheric  pressure  *  this  value,  used  for  transducer 
230  REM  pressure  calculations  on  any  transducer  which  is  not 

240  REM  referenced  to  vacuum,  can  be  read  off  any  RELIABLE 

250  REN  source  (i.e.  the  NKS  PVS-2  Portable  Vacuum  Standard, 

260  REM  or  the  Room  148  wall  barometer  if  it.  haa  been  calibrated). 

270  REN 

280  REN  4)  The  Vacuum  Reference  Pressure  -  this  Is  the  vacuum  reference 
290  REM  pressure  applied  to  the  reference  loort  of  transducers  placed 

300  REM  in  the  test  section  area.  It  is  r«tsd  off  the  20  me  Hg  range 

310  REM  vacuum  gauge  which  reflects  vacuum  j^ressure  at  the  tips 

320  REM  of  the  vscuue  reference  tubing  (not  the  vacuus  pump). 

330  REN 

340  REN  5)  Heliue  Line  2nd  Dome  Pressure  *  this  is  the  pressure  setting  in 
350  REM  paig  for  the  particular  run  of  the  2nd  Grove  dome  regulator 

360  REM  gauge  located  on  the  heliue  farm. 

370  REM 


380  REM  6)  Diffuser  Throat  Diameter  -  this  is  the  measurement  in  inches 
390  REM  of  the  supersonic  diffuser  throat  diameter. 

400  REM 

410  SCREEN  9:  COLOR  10,1:CLS  0:PRINT:PRINT 
420  PRINT  SPC(28)  "PRE-RUN  INPUTS  PROGRAM* 

430  PRINT  . . . 

440  COLOR  7,1  4 

460  COUNTERS  «  "D:\CB\DATAFILS\COUNTER.DAT* 

470  OPEN  COUNTERS  FOR  INPUT  AS  «1 

480  I NPUTS1, COUNTS 

490  CLOSE  i1 

500  COUNTS  >  COUNTS^I 
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S10  OPEN  COUNTERS  FOR  OUTPUT  AS  #1 
520  URITE#1,CQUNTX 
530  RUNNUMX  >  COUNTS 
540  CLOSE  #1 

550  PRINT  SPC(15)  "PLEASE  ENTER  THE  FOLLOWING  PRE-RUN  «";RUNNUMX;"  INPUTS:" 

560  PRINT 

570  RUNNUHSsSTRS(RUNNUMX> 

580  N  >  LEN(RUNNUM$) 

590  INPUT  "  HeliuM  Line  Temperature  in  degrees  F:  ",HETEHPI 

600  HETEMPI>HETEHPI+459.67  'Convert  to  degrees  R 

610  INPUT  "  Air  Line  Temperature  in  degrees  F:  ".AIRTEMPI 

620  AIRTEMP!>AIRTEMPIH59.67  'Convert  to  degrees  R 

630  INPUT  "  Atmospheric  Pressure  in  inches  of  Hg:  ",PATNI 

640  PATMI>  PATH)*  .49124  'Convert  to  psia 

650  INPUT  "  Vacuus  Reference  Pressure  in  .ns  Hg:  ",PREFI 

660  PREFI  >  PREF!  *  .01933  'Convert  to  psia 

670  INPUT  "  Helius  Line  2nd  Dook  Pressure  in  psig:  ",D0ME2PRESSI 

680  INPUT  "  Diffuser  Throat  Oiaswter  in  inches:  ",DIFFDIAI 

690  PRERUNS  »  "O:\CB\DATAFILS\PRE"  ♦  RIGHTS(RUNNUMS,N-1)  ♦  ".OAT" 

700  OPEN  PRERUNS  FOR  OUTPUT  AS  #1 

710  WRITE  »1,HETEMPi,AIRTENPI,PATNI,0aME2PRESSI.0IFF0IAI 
720  CLOSE  #1 

730  COLOR  10.1 :CLS  0:PRINT:PR1NT:PRINT:PRINT:PRINT:PRINT:PRINT: 

740  INPUT"  PRESS  CENTER]  TO  EXIT.  THEN  BEGIN  NICOLET  500  DATA  ACQUISITION  PROGRAM". DUMMYS 
750  COLOR  7.1:CLS  0 
760  SYSTEM 
770  END 


)8TRUW.BA8 


10  '  PROGRAM  POSTRUN. BAS 
20  ' 

30  '  Written  by:  Curt  D.  Botts 

40  '  This  program  reduces  wavefons  data  saved  in  Nicolet  .WFT  files,  iuo  output 

50  '  files  are  created  for  each  channel  for  each  run.  A  "(run  nu(ber)00(channel  nunber).DAT" 

60  '  file  is  compatible  with  the  GRAPHER  program  for  creating  pressure  versus  time  graphs.  A 
70  '  "(run  nupber)OUT.DAT"  file  contains  information  on  a  single  run  including  nunerical  output 
80  '  of  pressures,  mass  flows,  gaaam  value,  and  exit  Mach  number  at  time  intervals  during  a  run. 
90  ' 

100  OPTION  BASE  1 

110  DIM  SLOPEI(20).PRESSI(20.101).HOOTAIR((101).MDOTHELI(101).GAMMAI(101).HOLDI(20}.MACHEI(101) 
120  ON  ERROR  GOTO  140 
130  GOTO  240 
140  A«ERR:B«ERL 

150  SCREEN  0:COL0R  28.1 :CLS  0:PRINT:PRINT:PRINT 

160  PRINT  SPC(19)  "ERROR  NUMBER  ".-A; "OCCURRED  AT  LINE";B 

170  COLOR  7.1 

180  PLAY  ■03":PLAY"L4";PLAY  "C":PLAY"L6":PLAY  "C":PLAY"L8":PLAY  "C":PLAY"L4":PLAY  "C" 

190  PLAY"L4";PLAY"E-":PLAY"L6":PUY"0":PLAY"0":PUY"C";PUY"C":PLAY"L4";PLAY"02" 

200  PLAY"B";PLAY"Q3":PUY"L1":PUY"C" 

210  GOTO  2250 

220  • 

230  REM  Input  transducer  calibration  slopes. 

240  DATA  0.87693.0.592344.0.82331.2.49835.17.1952.0.924565.0.977901.0.700219 
250  DATA  0.931072.0.897628.0.782386.0.437291.0.786606.0.7800330.0.988232.2.56069 
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CLOSE  #1 


910 
920  ' 

930  REM  Retrieve  zero  voltage  values  froai  ZER001A.UFT  through  ZER005D.UFT 
940  REN  data  files  by  first  converting  theai  to  .FLT  files  via  the  WFT2FLT.EXE 
950  REM  prograia.  This  prograia  reads  the  .WFT  binary  data  and  creates  a 
960  REN  2  coluan  (volts, tiaie)  .FLT  file.  The  average  over  all  data  points 
970  REM  per  channel  becowes  that  channel's  zero  voltage. 

980  ZEROS^'H) :  \CBVDATAF  I  LS\ZER00"-»RI  GHTS(CHANNEL$,2) 

990  SHELL  ■n):\CB\DATAFILS\UFT2FLT  ZER(»*".WFT>JUNKOUT“ 

1000  OPEN  ZEROS^^.FLT"  FOR  INPUT  AS  #1 

1010  REN  This  finds  the  first  data  point. 

1020  FOR  JX>1  TO  20 

1030  INPUT  #1,CHECKI 

1040  IF  CHECKIo.06  THEN  HOLDI ( JXXHECKI 

1050  IF  CHECKI-.06  GOTO  1090 

1060  NEXT  JX 

1070  PRINT  "  ERROR  IN  READING  DATA  FILE  ZEROS; “.FLT" 

1080  GOTO  180 

1090  VOLTSUMI=HOLOI(JX-1)>HOLOI(JX-3) 

1100  FOR  JX>1  TO  998 

1110  INPUT  «1, VOLT I, SECS I 

1120  VOLTSUMI-VOLTSUMI+VOLTI 

1130  NEXT  JX 

1140  ZEROVOLTI>VOLTSUNI/1000l 

1150  CLOSEAI 

1160  • 

1170  REM  Convert  .WFT  file  to  .FLT  file  via  the  Nicolet  WFT2FLT.EXE  program. 

1180  REN  This  reads  the  .UFT  binary  data  and  creates  a  2  coltnn  (volts, secs) 

1190  REN  .FLT  file. 

1200  SHELL  “O:\CB\DATAFILS\UFT2FLT  DATAFILES>“>JUNKOUT” 

1210  OATAF I  LES>“D :  \CB\DATAF  I  LS\“+RUNNUMS4^"00“+RIGHTS(CHANNELS,  2) 

1220  REM  To  save  disk  space  write  over  old  .FLT  file. 

1230  NAME  0ATAFILES+".FLT“  AS  “O:\CB\OATAFlLS\FLTHOLD.OAT" 

1240  OPEN  “O:\CB\DATAFILS\FLTHOLD.DAT"  FOR  INPUT  AS  «1 

1250  REM  Open  .DAT  file  to  store  pressure  vs.  time  data  which  is  accessible 

1260  REM  from  GRAPHER. 

1270  OPEN  DATAFILES4^".DAT“  FOR  APPEND  AS  «2 

1280  REN  This  finds  the  first  data  point. 

1290  FOR  JX*1  TO  10 

1300  INPUT  «1, CHECK! 

1310  IF  CHECKIo.001  THEN  LASTI>CHECK[ 

1320  IF  CHECKI>.001  GOTO  1380 

1330  NEXT  JX 

1340  PRINT  "  ERROR  IN  READING  DATA  FILE  ”;DATAFILES;".DAT" 

1350  GOTO  2250 

1360  • 

1370  REM  Calculate  the  pressure  at  each  10th  time  point  and  store  in  array  PRESS(channel,time). 
1380  NX>0:CX-1 :NX>0 

1390  FOR  JX>1  TO  1000 

V,00  NX-NX^1 

1410  IF  JX>1  THEN  VOLTI-LASTI :SECS|sCHECKI :GOTO  1450 

1420  INPUT  «1, VOLT  I, SECS I 

1430  IF  JX>1000  GOTO  1450 

1440  IF  NXolO  GOTO  1540 

1450  NXsMX^I 

1460  DELVOLTI-ABS(VOLTI-ZEROVOLTI) 

1470  PRESS!  (IX,NX)>OELVOLT!*SLOPEI(IX)-'^PREFI 

1480  REN  Except  those  referenced  to  atmosphere  arxi  differentials  without  reference. 

1490  IF  (KX34)  AND  (LX>68)  THEN  PRESS! (IX,NX}«DELVOLT!*SLOPE! (IX) 

1500  IF  (KX>5)  AND  (LX-65)  THEN  PRESSI(IX,NX)«OELVOLTI*SLOPEI(IX) 

1510  IF  (KX>5)  AND  (LX-66)  THEN  PRESS! (IX,NX)>PRESS!(IX,MX)+PATN! 

1520  NX>0 

1530  WRlTEiit2,SECS!, PRESS!  (IX,NX) 

1540  NEXT  JX 

1550  CLOSE  «1:CL0SE  «2:SHELL  "DEL  D:\CB\0ATAFILS\FLTH0L0.DAT” 

1560  COLOR  2,1; PR I NT"  COMPLETED": COLOR  7,1 

1570  NEXT  LX 
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'Calculate  the  air  naas  flow  rate. 
■Calculate  the  heliui  Maaa  flow  rate 


1S60  NEXT  ICX 

1590  • 

1600  GOSUB  2260 

1610  G0SU8  2420 

1620  ■ 

1630  REN  Calculate  the  Mixture  gaeim  value  and  Che  nozzle  exit  Nach  nuaber. 

1640  REN  Write  to  GRAPHER  coMpatible  files. 

1650  OPEN  ••d:\cb\datafil8\">RUNNlJ*a+‘'gaaMW.clat''  FOR  APPEND  AS  «1 

1660  OPEN  •'d:\cb\datafils\''-MlUNNlJMS««sache.dat'‘  FOR  APPEND  AS  «2 

1670  FOR  CX>1  TO  100 

1680 

GANNAI (CX}*(NDOTAIRI (CX)*CPAIRIiNDOTHELI (CX)*CPHELI )/(ND0TAIRi (CX)*CVAIRI«NDOTHELI (CX)*CVHELI ) 

1690  NACHEI  (CX)>saR(2l/(GMMAI  (CX)-1 1  )*«(PRESSI(5,CX>/PRESSI (15,CX))*(1 1  -1  l/GANNAI  (CX)))-1 1 )) 

1700  PRINT#1,GANNAI(CX),CSNG(CX-1)/10i 

1710  PRINT#2,NACHEI(CX),CSNG(CX-1)/10I 

1720  NEXT  CX 

1730  CLOS£#1:CLOSE#2 

1740  ' 

1750  REN  Create  output  data  file  for  this  run. 

1760  OPEN  ■A'',#1,'n):\C8\DATAFaS\“+RUNNUN**'WT.0AT" 

1770  PRINT#1,SPC<28)  "RUN  NIMBER  ••;RUNNUtt;“  0ATA“:PRINT#1,'"';PRINT#1,“" 

1780  PRINTRI,"  Date:  ■; 

1790  IF  NONTH<10  THEN  PRINTRI, USING  "0«  /".-NONTH; :GOTO  1810 
1800  PRINTRI, USING  "«#  /"; NONTH GOTO  1810 

1810  IF  0AY<10  THEN  PrTnTRI .USING  "O#  /“.-DAY; :GOTO  1830 
1820  PRINT#1, USING  "R#  /".'DAY; 

1830  PRINTRI, USING  “RR^.'YEAR 

1840  HOURS«SECONOS/3600000I:IHQURS>INT(HOURS} 

1850  NINUTES-(HOURS*IH0URS)*60:ININUTES>INT(NINUTES) 

1860  SEC0N0S>(NINUTES-ININUTES)*60 
1870  PRINTR1  ”  Tine:  ■** 

1880  IF  IHOURS<10  THEN  PRINTRI  .USING  "O#  I  HOURS; :  GOTO  1900 
1890  PRINTRI. USING  "RR  :'';IHOURS; 

1900  IF  ININUTES<10  THEN  PRINTRI .USING  "OR  :";ININUTES;:GOTO  1920 
1910  PRINTRI. USING  "RR  :";ININUTES; 

1920  IF  SECON0S<10  THEN  PR1NTR1.USINC  -OR.«R";S£CONDS:GOTO  1940 

1930  PRINTRI. USING  "RR.RR";SECONOS 

1940  PRINTRI. USING  "  Atnospheric  pressure.:  M.RRR  psia";PATNI 

1950  PRINTRI. USING  "  Reference  vacuus  pressure  :  RR.RRR  psia";PREFI 

1960  PRINTR1.USING  •  Nelius  tenperature  :  m.m  deg  F";HETENP I -459.67 

1970  PRINTRI. USING  "  Air  tenperature  :  m.m  deg  F";AIRTENPI -459.67 

1980  PRINTRI. USING  ■  Heliua  2nd  done  pressure  :  RRR.RR  p6ig";D0ME2PRESSI 

1990  PRINTRI. USING  "  Diffuser  throat  dimeter":  RR.RRRR  inch";DIFFDIAI 

2000  PRINTRI."" 

2010  FOR  JX>1  TO  100 

2020  PRINTRI," _ 

2030  PRINTRI, USING"  Tine«R.RR  sec";CSNG(JX-1)/10l:PRINTR1,"" 

2040  PRINTRI,"  psia  at  positior#:” 

2050  PRINTRI,"  . « 

2060  PRINTRI, USING"  1)  RRR.RRR  12)  RRR.ihINf  ";PRESSI(1,JX},PRESSI(9,JX} 

2070  PRINTRI, USING"  2)  RRR.RRR  13}  RRR.RR#  ";PRESSI (2, JX),PRESSI (10, JX) 

2080  PRINTRI, USING"  3)  RRR.RRR  14)  RRR.RRR  Air  nass  flow  rate^RR.RRR 

lbm/s";PRESSI(3,JX),PRESSI(11,JX),MDOTAIRI(JX) 

2090  PRINTRI, USING"  5)  RRR.RRR  15}  RRR.RRR  Heliua  nass  flow  rste-RR.RRR 

lbm/s";PflESSI (4, JX),PRESSI (12, JX),MDOTHELI ( JX) 

2100  PRINTRI, USIi«G"  6)  RRR.RRR  16)  RRR.iMI  Mixture 

ganma>RR.RRR”,-PRESSI(5,JX),PRESSt(13,JX),GM«M(JX) 

2110  PRINTRI, USING"  9)  RRR.RRR  17)  RRR.RRR  •;PRESSI(6,JX).PRESSI(14,JX) 

2120  PRINTRI .USING"  10)  RRR.RRR  18)  RRR.RRR  Nozzle  exit 

Msch>RR. RR" ; PRESS I ( 7, JX) , PRESS I ( 1 5 , JX) , MACHE I ( JX) 

2130  PRINTRI, USING"  11)  RRR.RRR  31)  RRR.RRR  ";PRESSI (8, JX),PRESSI (18, JX) 

2140  NEXT  JX 
2150  CLOSE  R1 

2160  PUY"M8":PLAY"L4":PLAY"03";PLAY"C":PLAY"0";PUY"C":PUY"0";PLAY"C":PLAY"0" 

2170  IF  CHOICEX*1  GOTO  2190  ELSE  NEXT  RUNNUMX 
2180  GOTO  2230 

2190  INPUT"  WOULD  YOU  LIKE  TO  REDUCE  ANOTHER  RUH";ANSUER$ 
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2200  IF  ANSUER^"YES"  OR  ANSUERS>>”y«*”  OR  ANSUERR-^V*  OR  ANSWERS-^y*  GOTO  460 
2210  IF  ANSUERS-^NO*  OR  ANSWERS-^no*  OR  ANSUERS^^N"  OR  AkSUBRU^n"  GOTO  2230 
2220  PRINT"  PLEASE  ANSWER  YES  OR  NO":GOTO  2190 

2230  COLOR  7.1 
2240  SYSTEM 
2250  END 

2260  REM  Mass  flow  rat*  calculation  through  air  orific*  aeter. 

2270  RI>1717.S9  'air  gaa  constant  (sq.ft/sq.sec-deg  R) 

2280  Kl«1.4  'air  ratio  of  specific  heats 

2290  BETA) >.2666667  'ratio  of  orifice  to  pipe  diaaieters 

2300  FOR  MX>1  TO  100 

2310  IF  PRESSI(17,MX)<>0I  GOTO  2330 

2320  PRINT:PRINT  SPCdS)  "BAD  RUN!  FLOW  DID  NOT  OCCUR  AT  TRIGGER  I  ":RETURN 
2330  RATI01>(PRESSI(18,NX)-PRESSI(17,NX))/PRESSI(18,MX) 

2340  IF  RATIO-11  THEN  MDOTAIRI (NX)-OI :GOTO  2400 

2350  PART1l»RATI0l*(2l/KI)*KI/(ICI-1l)"(1l-(RATI0r((ICI-1l>/ICI>)>/(1l'RATI0l) 
2360  PART2l-PART1l*(l!-B£TAI*4l)/<1l-BETAr4l*RATI0r{2l/KI)) 

2370  YI«SQR<PART2i) 

2380  PARTI l-SQR(PRESSI (18, MX)/(RI*AIRTEMPI)*PRESSI(17.NX)) 

2390  MDOTAIRI (NX}-13.760371«"YI*PART1 I 
2400  NEXT  NX 
2410  RETURN 

2420  REN  Mass  flow  calculation  through  Helius  venturi  aster. 

2430  RR- 12427. 9291#  'helius  gas  constant  (sq.ft/sq.sec-deg  R) 

2440  Xl-I .666667  'helius  ratio  of  specific  heats 

2450  BETAI-. 59809  'rstio  of  throat  to  pipe  diaaeters 

2460  FOR  NX-1  TO  100 
2470  IF  PRESSI(16,NX)<>0I  GOTO  2430 

2480  PRINT:PRINT  SPC(18}  "BAD  RUNI  FLOW  DID  MOT  OCCUR  AT  TRIGGERI":RETURN 
2490  RATI0I»(PRESSI(19,NX)-PRESSI(16,NX))/PRESSI(19.MX) 

2500  IF  RATIO-11  THEN  MDOTHELI (MX)-OI :GOTO  2570 
2510  IF  DOHE2PRESSI-0I  THEN  MDOTHELI (MX)-O I : GOTO  2570 

2520  PART1l-RATI0r(2l/KI)*KI/(KI-1l)*<1l-(RATI0r((KI-1l)/KI)))/(1l-RATI0l) 
2530  PART2l-PART1l*(1l-BETAr4l)/<1l-BETAr4l*RATI0r(2l/KI)> 

2540  YI-S0R<PART2l ) 

2550  PART 1 1 -SOR ( ( PRESS I ( 19, MX)/ ( R I "HETEMP I ) )*PRESS I ( 16, MX) ) 

2560  MDOTHELI <MX)*9.5272834#*YI*PART1 I 
2570  NEXT  MX 
2580  RETURN 
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Mixing  of  primary  cold  flow  air  and  secondary  helium  to  control  the  ratio  of 
specific  heats  for  the  medium  flowing  through  a  simulated  chemical  laser 
nozzle/lasing  cavity  was  accomplished.  The  effects  of  a  range  of  mixture 
specific  heat  ratios  on  flowfleld  behavior  were  examined  using  static  pressure 
ports  in  the  lest  cavity.  Schlieren  photography  and  high  speed  filming  aided 
description  of  the  flow  dynamics.  Results  indicated  that  boundary  layer  effects 
became  evident  in  the  nozzles  as  specific  heat  ratios  increased.  Large  pressure 
fluctuations  were  observed  in  the  cavity  when  helium  was  introduced  into  the  flow 
to  raise  the  specific  heat  ratio.  This  unstable  behavior  was  attributed  to  the 
helium  mass  flow  into  the  mixer  and  the  mixer  design  Itself.  Use  of  the 
air/helium  mixer  brought  about  the  pressure  fluctuations  earlier  in  a  test  run 
than  with  the  mixer  removed  under  the  same  conditions.  Favorable  pressure 
conditions  for  lasing  were  achieved  for  at  least  two  seconds  for  the  supersonic 
nozzles’  design  specific  heat  ratio  of  1.51.  Adverse  pressure  behavior  was  also 
attributed  to  three  dimensional  viscou.;  effects  along  the  cavity  walls. 
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